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1.  INTRODUCTION 

The  1980  Monterey  Bay  Turbulence  model  verification  experiment 
was  conducted  to  validate  marine  surface  layer  turbulence  models 
to  be  used  with  the  Navy  High  Energy  Laser  (HEL)  program.  The 
experiment  consisted  of  laser  scintillometer  measurements  across 
Monterey  Bay  accompanied  by  simultaneous  shipborne  bulk  meteorology 
and  turbulence  measurements.  In  addition,  an  aircraft  was  flown 
to  provide  turbulence  and  meteorological  measurements  from  near 
the  surface  to  above  the  marine  inversion,  and  GOES  and  DMSP 
satellite  imagery  was  obtained  in  the  visible  and  infrared  channels. 

The  experiment  verified  that  the  marine  surface  layer  turbu¬ 
lence  model  recommended  to  the  HEL  program  by  NEPRF  (Burk,  Goroch, 
Weinstein  and  Panofsky,  1979)  agreed  with  the  optical  measurements 
to  a  factor  of  two.  Experimental  results  indicate  that  the  model 
depends  on  conventional  bulk  sea  surface  temperatures  (platinum  or 
bucket  measurements)  rather  than  infrared  temperatures.  Use  of 
infrared  temperatures  in  the  model  did  not  show  any  agreement  with 
optical  measurement.  It  is  possible  that  a  bias  might  be  added  to 
IR  measurements  to  provide  better  agreement. 

2.  BACKGROUND 

Turbulence  in  the  atmosphere  affects  optical  signal  propaga¬ 
tion  by  changing  the  resulting  beam  characteristics.  The  effect 
of  turbulence  is  characterized  by  the  optical  index  of  refraction 
structure  function,  Cj|,  which  is  related  to  ambient  meteorologi¬ 
cal  conditions.  The  HEL  program  requested  NEPRF  to  conduct  a 
workshop  to  recommend  a  model  correlating  with  meteorological 
parameters . 

The  recommended  model  is  discussed  thoroughly  in  the  report 
of  the  workshop  (Burk,  et_  a_l_.  ,  1  979).  Briefly  summarized,  the  model 
modifies  the  bucket  sea  surface  temperature  for  surface  effects 
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T(bucket)  -  0.025  T  .  -  0.1  (V  <  6  m  sec"  ) 

o  o  all 
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T  (bucket)  -  0.1  T  .  -  0.8  +  V ( 0 . 1 1 7  + 0 . 01 25  T  .  ) 
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where  Tss  (bucket)  is  the  bucket  sea  surface  temperature,  T  . 
is  air  temperature,  and  V  is  10  m  wind  speed.  'The  following 
formulas  provide  the  surface  layer  momentum  heat  and  moisture 
fluxes  in  the  form  of  scaling  wind  U*  temperature;  T*  and  specific 
humidity  Q  * 

U*  =  kU[ln(z/z0)  -  ^n(z/L)]_1 

Q*  =  |  (q  -  qs)[ln(z/zQq)  -  ( z/ L ) ] -1 

T*  =  |  (T  -  Ts)[ln(z/zoT)  -  ^(z/L)]"1 

zQ  -  "roughness"  parameter  for  momentum  2.  0.6  mm 

zQp  -  "roughness"  parameter  for  humidity 

zqj  -  "roughness"  parameter  for  temperature 

R  -  0.74 


-  von  Karmann's  constant  0.35 
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The  roughness  parameters  describe  the  detailed  molecular 
interaction  between  the  sea  surface  and  the  adjacent  air.  Their 
calculation  is  described  in  Burk  et  al . ,  1980. 

The  equations  provided  above  are  dependent  on  the  Monin- 
Obukh'ov  length  L.  The  ratio  of  height  z ,  to  L  is  given  by 


z_ 
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gkT* 

- 2(1 

(T  •  +  273 . 1 6 ) U* 


0.07 
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where 
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o 
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acceleration  due  to  gravity  980  cm  sec 

c  T* 

Bowen  ratio  = 

lH  w* 

specific  heat  of  air 

latent  heat  of  evaporation  of  water 


The  above  equations  are  solved  by  iteration.  An  initial 

3  0 

value  is  assumed  for  L  (L  =  1  •  1 0m  ).  The  scaling  parameters  U*, 
T* ,  q*  are  calculated  and  L  is  recalculated  and  compared  to  the 
previous  estimate.  When  two  estimates  agree  to  better  than  0.1%, 
the  temperature  structure  function  parameter  is 
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2 
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g(z/L) 


4.9(1  +2.75  y)  L  >_  0 
4.9(1  -  7  )“2/ 3  L  <  .0 


and  the  optical  index  of  refraction  parameter  is 


'T 


7 .81 7 • 1 0~5  p 

(T  .  +  273.1  5) 

v  air  ' 


r(l  + 


.03 
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where  p  is  air  pressure  in  millibars 


This  model  was  used  with  the  surface  correction  and  without 
for  bulk  meteorological  variables,  with  platinum  wire  and  infrared 
sea  surface  temperature  measurements  alternately. 
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3.  EXPERIMENTAL  PROCEDURE 

The  experiment  consisted  of  two  types  of  measurements:  shore 
to  shore  optical  scintillation  measurements,  and  shipborne 
meteorology  measurements.  The  measurement  techniques  are  out¬ 
lined  below;  complete  descriptions  are  included  in  the  reports 
of  Crittenden,  Milne,  Rodeback  and  Kalmback  (1980)  and  Davidson, 
Schacher,  Fairall,  Crittenden  and  Milne  (1980). 

4.  OPTICAL  MEASUREMENTS 

The  optical  equipment  consisted  of  a  DF  laser  operating  at 
10.6  micrometers  wavelength  situated  on  the  eastern  shore  of 
Monterey  Bay.  A  chopped  signal  was  sent  along  a  13.16  km  path 
over  Monterey  Bay  to  a  receiving  station  on  the  south  west  shore¬ 
line.  The  average  height  above  the  sea  surface  was  14  m.  The 
optical  range  is  shown  in  Figure  1. 

The  signal  was  measured  for  one  minute  to  provide  a  distribu¬ 
tion  of  signal  intensity.  The  standard  deviation  of  the  distribu- 

2  2 
tion  cfj  is  related  by  C^  by 

C2  =  1  .041  10"14  a  T 

n  I 

The  measurement  is  susceptible  to  saturation,  that  is,  an 

increase  of  turbulent  intensity  will  not  result  in  an  increased 

standard  deviation  of  the  intensity  distribution.  The  minimum 

intensity  value  for  the  occurrence  of  turbulence  has  been  suggested 

2 

by  Crittenden  to  be  at  o  =  0.8.  For  this  value,  the  corresponding 

2  - 1  5 

turbulent  intensity  is  Cn  =  6.6-10  This  value  was  not  reached 

during  the  experiment. 

5.  SHIPBORNE  MEASUREMENTS 

The  shipborne  measurements  consisted  of  bulk  meteorological 
parameters  and  derived  turbulence  parameters.  The  latter  measure¬ 
ments  were  taken  to  provide  an  indication  of  validity  of  the  thermal 

2 

and  humidity  components  of  the  bulk  model,  if  the  model  C  did 
not  agree  with  optical  measurements. 
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Figure  1.  Map  of  Monterey  Bay  and  optical  range. 


5 


The  bulk  measurements  consisted  of  cup  anemometers  and  wind 
vanes  for  wind  speed  and  direction.  These  were  corrected  for 
ship  motion.  Air  temperature  was  measured  at  4.2  m,  7.0  m  and 
19.6  m  above  the  surface.  Dew  point  was  measured  at  7.0  m,  and 
humidity  was  measured  at  19.6  m.  Sea  surface  temperature  was 
measured  by  a  platinum  thermometer  maintained  in  the  first  25  cm 
below  the  sea  surface.  The  infrared  sea  surface  temperature  was 
measured  by  an  infrared  radiometer  mounted  on  a  3  m  boom  on  the 
ship’s  rail. 

6.  SATELLITE  IMAGERY 

GOES  satellite  imagery  was  stored  on  magnetic  tape  for  the 
dates  28  April  through  9  May  1980.  The  (western)  GOES  data  was 
received  at  NEPRF  at  1130  (PDT)  daily  and  transferred  to  tape. 

The  data  was  in  the  form  of  visual  and  infrared  (8-12  pm)  images, 
with  256  gray  shades,  and  a  resolution  of  .3  n  mi. 

The  data  was  later  analyzed  as  follows.  Visual  imagery  was 
screened  to  determine  which  days  were  cloud  free  in  the  Monterey 
Bay  vicinity.  The  corresponding  IR  data  was  magnified  by  a  factor 
of  4  and  displayed  on  the  NEPRF  SPADS  (Satellite  Processing  and 
Display  System).  The  gray  shade  contrast  settings  were  adjusted 
until  the  sea  surface  temperature  contrast  was  maximized.  Two 
calibrations  were  used  to  relate  temperature  to  gray  shade  -  the 
first  was  a  general  calibration  accurate  to  1.0  degrees  C,  the 
second  accurate  to  0.5°C.  The  two  resulting  analyses  are  shown 
in  Figures  2  and  3  respectively. 

DMSP  data  was  also  collected.  Several  examples  of  anomalous 
gray  shades  occurred,  where  the  DMSP  image  (0.4  -  1.1  pm  sensitive 
range)  shows  contrasting  regions,  without  cloud  cover.  The  explana¬ 
tion  of  this  phenomenon  has  been  attributed  to  dryer  regions, 
presence  of  smaller  (drier)  aerosols  and  different  surface  ref¬ 
lectivity  (Fett  and  Isaacs,  1979). 

A  DMSP  image,  taken  on  1  May  1980,  shows  a  pronounced  darker 
shade  in  the  northern  part  of  Monterey  Bay.  The  darker  shade 
indicates  that  less  signal  is  being  received  at  the  detector.  The 
shape  of  the  darker  region  suggests  that  the  shade  is  wind  related. 
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Figure  3.  Analysis  of  GOES  satellite  data  7  May  1980. 


7 


Several  possible  mechanisms  attenuate  the  signal.  These  include 
decrease  of  large  aerosol  concentration  by  drying  overland  and 
smoothing  of  the  sea  surface  by  the  barrier  of  land  upwind  of 
the  darker  portion.  The  image  is  included  as  Fig.  4.  These 
phenomena  are  being  studied  more  extensively  to  explain  the 
uncertainties  of  the  image  (Fett,  R.,  personal  communication). 

7.  WEATHER  CONDITIONS 

The  basic  weather  conditions  occurring  during  the  entire 
experimental  period  were  neutral.  The  neutral  conditions  were 
modified  by  several  occurrences  of  moderately  stable  and  moderately 
unstable  conditions. 

The  upper  air  pattern  is  shown  on  the  National  Meteorologi¬ 
cal  Center  500  mb  charts  (Figures  5-17).  The  California  coastal 
region  was  characterized  by  a  weak  low  pressure  area  remaining 
stationary  from  before  the  experiment  to  6  May  1980.  This  stationary 
low  was  circled  by  serveral  small  weak  cold  fronts.  Thus  the 
synoptic  flow  was  in  general  weak,  with  alternating  periods  of 
onshore  and  offshore  flow.  The  onshore  flow  was  generally  accom¬ 
panied  by  low  level  clouds,  while  offshore  flow  was  cloud  free. 

These  conditions  started  to  be  modified  by  the  presence  of 
a  strong  low  pressure  area  off  the  British  Columbia  coast.  This 
was  responsible  for  the  clearer  conditions  at  the  end  of  the 
experiment.  The  cold  front  passed  the  Monterey  Bay  area  after 
the  experimental  period  ended. 

8.  RESULTS 

o 

The  index  of  refraction  structure  function,  C  ,  measured  by 

2  ^ 

the  optical  means,  was  compared  to  the  C  predicted  by  the  workshop 

2  ^ 

model.  The  range  of  optical  Cn  was  limited  to  relatively  low 
values  because  of  the  neutral  and  near  neutral  conditions. 

The  analysis  separated  the  data  into  three  groups,  as  sug¬ 
gested  by  the  Davidson,  et_  aj_  (1  980)  report.  These  groups  were 
(a)  unstable,  homogeneous  sea  surface  temperature  conditions,  (b) 
stable  homogeneous  sea  surface  temperature,  and  (c)  stable  condi¬ 
tions.  The  inhomogeneous  sea  surface  temperature  conditions  were 
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SUNDAY.  APRIL  27,  1W0 


Figure  5.  National  Meteorological  Center  surface  and  500  mb 
analyses,  27  April  1980, 
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MONDAY,  APRIL®,  1980 


HIGH  u 
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Figure  6.  National  Meteorological  Center  surface  and  500  mb 
analyses,  28  Apri 1  1980. 
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TUESDAY,  APRIL  29.  1980 


Figure  7.  National  Meteorological  Center  surface  and  500  mb 
analyses,  29  April  1980. 


WEDNESDAY.  APRIL  30. 1900 


Figure  8_  National  Meteorological  Center  surface  and  500  mb 
analyses',  30  April  1  980. 
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THURSDAY,  MAY  1,  1900 
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;re  9.  National  Meteorological  Center  surface  and  500  mb 
analyses,  1  May  1980. 
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FRIDAY,  MAY  2,  1980 


Figure  10.  National  Meteorological  Center  surface  and  500  mb 
anal yses  ,  2  May  1  980 . 


15 


SATURDAY,  MAY  3,  1W0 


Figure  11.  National  Meteorological  Center  surface  and  500  mb 
analyses ,  3  May  1 980 . 
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SUNDAY.  MAY  4,  1980 
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Figure  i  ^ .  National  Meteorological  Center  surface  and  500  mb 
analyses ,  4  May  1 980. 
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Figure  13.  National  Meteorological  Center  surface  and  500  mb 
analyses,  5  May  1 980 . 
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TUESDAY,  MAY  6,  1980 


Figure  14.  National  Meteorological  Center  surface  and  500  mb 
analyses,  6  May  1980. 
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WEDNESDAY,  MAY  7,  I960 


Figure  15.  National  Meteorological  Center  surface  and  500  mb 
analyses,  7  May  1 980 . 
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THURSDAY,  MAY  8,  1980 
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Figure  16.  National  Meteorological  Center  surface  and  500  mb 
anal yses ,  8  May  1 980. 
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Figure  17. 


National  Meteorological  Center  surface  and  500  mb 
analyses.  9  May  1980. 
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those  in  which  there  seemed  to  be  a  1 °C  to  2°C  warming  of  the 

sea  surface  temperature  under  the  optical  path  compared  to  the 

temperature  farther  from  shore  in  the  bay. 

Both  optical  and  model  results  had  a  certain  amount  of  un- 

2 

certainty  involved  in  the  measurements.  The  optical  Cn  had  an 
uncertainty  in  evaluating  the  standard  deviation  of  the  inten¬ 
sity  distribution.  This  uncertainty  showed  up  in  near  simul¬ 
taneous  measurements  having  a  maximum  of  about  15%  difference 
2 

in  .  Similarly  the  model  result  also  was  uncertain  because 

of  sea  surface  temperature  uncertainty.  The  uncertainty  in  air 

sea  temperature  difference  of  0.2  degrees  results  in  a  12% 

2 

uncertainty  in  the  value  of  C  . 

A  geometric  method  was  chosen  to  compare  logarithms  of  the 

2 

measured  and  model  values  of  C  .  The  method  basically  calculates 

the  perpendicular  "distance"  from  each  observation  to  the  line 

o 

corresponding  to  identical  model  and  measured  Cn  values  (see 
Fig.  18).  This  distance  was  calculated  for  all  measurements 
and  the  mean  d  and  standard  deviation  s  of  the  logarithm  of  the 
distribution  was  evaluated.  The  latter  two  numbers  provided  a 
quantitative  measure  of  how  well  the  optical  measurements  and 
the  model  provided  the  same  results. 

The  optical  data  were  compared  to  the  model  results  using 
platinum  wire  sea  surface  temperatures  ("bucket")  and  the  infrared 
derived  temperatures.  These  comparisons  used  sea  surface  temper¬ 
atures  with  and  without  the  sea  surface  corrections  described  in 
Section  2,  The  results  are  shown  in  Figures  19-22.  The  optical 
Cn  were  also  compared  to  Cn  evaluated  from  shipborne  Cj  measure¬ 
ments  in  Figure  23. 

9.  CONCLUSIONS 

The  comparison  of  optical  and  model  results  in  Figure  23  shows 

2 

good  agreement  of  the  two  sets  of  Cn  values.  In  particular,  the 
figures  indicate  that  66%  of  data  had  the  optical  data  within  a 
factor  of  2.1  of  the  model.  The  97%  and  99%  confidence  intervals 
were  factors  of  4.8  and  10.5  respectively. 

The  addition  of  the  sea  surface  temperature  correction  makes 
only  a  slight  difference  on  the  agreement  of  model  and  measurement. 
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d=R  *  sin(T— TT/4) 

R=l (x— x,)2  +  |y-y,)2  p 
T=arc  tan((y— y1)/(x— xl )) 

Figure  18.  Geometry  of  data  comparison. 
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Figure  20.  Comparison  of  optical  Cp  to  workshop  model  (uncorrected 
Pt  temperature) . 
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Figure  21.  Comparison  of  optical  to  workshop  model  (corrected 
I R  temperature). 
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Figure  22.  Comparison  of  optical  to  workshop  model  (uncorrected 
I R  temperature ) . 
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Figure  23.  Comparison  of  optical  Cn  to  Acania  turbulence  measurements. 


The  mean  and  standard  of  the  model  and  measured  differences 

changed  only  in  the  second  and  third  significant  figures.  Based 

on  these  results,  the  addition  of  the  temperature  correction  does 

? 

not  benefit  the  bulk  C  model. 

n  2 

The  comparison  of  optical  Cp  and  that  derived  from  infrared 

sea  surface  temperatures  shows  little  agreement  between  the  two. 

This  occurred  regardless  of  whether  the  surface  temperature 

correction  was  included.  The  compared  to  the  bucket  temperature, 

IR  measurement  seemed  to  have  a  bias  which  could  not  be  handled 

p 

by  the  model,  resulting  in  which  had  no  correlation  to  the 

optical  measurements.  Further  work  should  be  conducted  to  analyze 

the  reasons  for  the  bias  in  temperatures,  and  the  possibility 

of  formulating  a  model  of  sea  surface  temperature  bias  for  the 

bul k  C2  model  . 
n 

The  satellite  derived  IR  temperatures  were  in  agreement  with 
the  ship  IR  surface  temperatures,  more  than  with  the  bucket  tem¬ 
peratures.  The  development  of  a  bias  model  would  possibly  enable 

use  of  satellite  infrared  temperature  data  as  an  aid  to  mesoscale 
2 

Cn  evaluation. 

All  the  above  conclusions  are  based  on  this  data  set,  which 

was  limited  to  near-neutral  stability  conditions.  The  credibility 
2 

of  the  Cn  model  (based  on  land  measurements  and  model  development) 
is  best  in  unstable  conditions,  less  good  in  neutral  and  least 
in  stable  conditions.  Based  on  the  experimental  comparisons, 
the  bulk  model  is  shown  calculating  optical  with  accuracy  in 
neutral  conditions.  By  inference,  we  can  say  that  in  very  un¬ 
stable  conditions  the  model  calculation  would  be  better.  For  ulti¬ 
mate  verification  however  the  model  should  be  validated  in  very 
turbulent  conditions,  both  stable  and  unstable. 
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